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Abstract. Many X-ray accreting pulsars have a soft excess below 10 keV. This feature has been detected also in faint sources 
and at low luminosity levels, suggesting that it is an ubiquitous phenomenon. In the case of the high luminosity pulsars (Lx > 
10 36 erg s _1 ), the fit of this component with thermal emission models usually provides low temperatures (kT < 0.5 keV) and 
large emission regions (R > a few hundred km); for this reason, it is referred to as a 'soft' excess. On the other hand, we 
recently found that in persistent, low-luminosity (Lx ~ 10 34 erg s~ x ) and long-period (P > 100 s) Be accreting pulsars the 
observed excess can be modeled with a rather hot (/cTbb > 1 keV) blackbody component of small area (Rbb < 0.5 km), which 
can be interpreted as emission from the NS polar caps. In this paper we present the results of a recent XMM-Newton observation 
of the Galactic Be pulsar RX J0440.9+443 1 , which is a poorly studied member of this class of sources. We have found a best-fit 
period P = 204.96 ± 0.02 s, which implies an average pulsar spin-down during the last 13 years, with P ~ 6 x 10~ 9 s s _1 . 
The estimated source luminosity is Lx ~ 8 x 10 34 erg s" 1 : this value is higher by a factor < 10 compared to those obtained 
in the first source observations, but almost two orders of magnitude lower than those measured during a few outbursts detected 
in the latest years. The source spectrum can be described with a power law plus blackbody model, with kTBB ~ 1-34 ± 0.04 
keV and -Rbb = 273 ± 16 m, suggesting a polar-cap origin of this component. Our results support the classification of RX 
J0440. 9+4431 as a persistent Be/NS pulsar, and confirm that the hot blackbody spectral component is a common property of 
this class of sources. 
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1. Introduction 

Most of the X-ray binary pulsars (XBPs) are High Mass X- 
Ray Binaries (HMXRBs) in which a Neutron Star (NS) with 
magnetic field B ~ 10 12 G is accreting matter from a high- 
mass early-type star, either an OB supergiant or a Be star. They 
can be persistently bright, with luminosities in excess of 10 34 
erg s~ 4 , or transient sources characterized by quiescent phases, 
with emission around 10 34 erg s _1 or less, in terrupted by bright 
outbu rsts reaching Lx ~ io je> - ci8 erg s - 1 dNegueruelalli998t 
l lRridl2007tlsidol«2010h . 

In these sources the X-ray spectra between 0.1 and 10 
keV are usually described by a rather flat power-law, with 
photon index ~ 1, but several XBPs have shown a marked 
'soft ' X-ray excess above the main power-law component 
(see La Palombara & Mereghetti 2006 for a review); it is well 
described by a thermal emission model (either blackbody, 
bremsstrahlung or mekal) with low temperature (kT$E < 0.5 
keV) and large emission area (i?sE > a few hundred km). 
This feature has been detected not only in the high-luminosity 
sources (with Lx ~ 10 37-38 erg s -1 ) but also in several 



low-luminosity (Lx ~ 10 35-36 erg s _1 ) XBPs observed in 
the Small Magellanic Cloud (SMC), whe re its detection is 
favoured by the low interstellar absorption (Sasaki et al. 2003 
l Ueno et all l2004t iMaiid etail |2004 iHaberl & Pietschl £005 



lHaberl et al.ll2008 ). Only in a few cases this low-energy com- 
ponent showe d coherent pulses an d the debate over its origin 



remains open. Hickox et al. (2004) have shown that a soft spec- 
tral component could be a very common, if not ubiquitous, fea- 
ture intrinsic to X-ray pulsars: it is visible in all sources with 
a sufficiently high flux and small absorption, and its origin is 
related to the source total luminosity. 

Recently, based on XMM-Newton data, we have 
observed a clear thermal excess also in three of the 
four persiste n t Be pulsars originally identified by 
teig&Rochd dl999t). i.e. RX .T0146.9+6121/LS I +61° 235 

4U0352+309/XPersei 
20071) 



2006), 



and 



dLa Pa lombara et alj l2009h . 



(La Palombara & Mereshett 
( La Palombara & Mereghett 
RX J1037.5-5647/LS 1698 
These three sources are characterized by a persistently low 
luminosity (Lx ~ 10 34 ~ 35 erg s _1 ) and a long pulse period 
(P > 100 s). These properties suggest that the NS orbits 
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the Be star in a wide and nearly circular orbit, continuously 
accreting material from the low-density outer regions of 
the circumstellar envelope; in the case of 4U 0352+309 this 
picture is supported by the long orbital period of 250.3 days 
(Del gado-Martf et al.l 1200 ll) . For these sources the detection 
of the thermal component was favoured by the small distance 
(d < 5 kpc) and interstellar absorption (TVh ~ 10 21 cm -2 ). We 
found that their soft excess can be fitted only with a blackbody 
(other simple models are rejected), which contributes for 
30-40 % of the total flux; interestingly, in comparison with the 
other, more luminous sources, their blackbody component is 
characterized by a higher temperature (/cTbb > 1 keV) and a 
much smaller emission radius (Rbb < 0.5 km). This hot BB 
spectral component sets these low-luminosity and long-period 
sources apart from all the other pulsars, strongly suggesting 
that they form a distinct class. In their case the thermal 
component could be due to a different emission mechanism 
than in the high-lu minosity pulsars. Based on the work of 
Hic kox et al.1 (120041) . it can be attributed to emission from the 
neutron-star polar caps. This is supported by the emission 
area of the blackbody component, which is consistent with the 
estimated polar cap size, and by the fact that the low energy 
part of the spectrum is clearly pulsed. 

In this paper we present the results of a recent XMM- 
Newton observation of RX J0440. 9+443 1 , the remaining mem- 
ber of this class of Be/NS pulsars. This system was discov- 



Marc h 2010 dMorii et alj|2010l) and b y RossiXTE on 6 April 
2010 dFinger & Camero-Arranzl2010l) ; the seco nd one was ob- 



ered during the ROSAT Galactic plane survey (Mo tch et al 



19971) and identified with LS V +44 17, a moderately red- 
dened (E(B - V) = 0.65 ± 0.05) B0.2 Ve star at - 3.3 kpc 
(Reig 2011). Thanks to observations with the PCA instrument 
on board RossiXTE. iReig & Roche (1999) performed the first 
detailed timing and spectral analysis, and discovered a pulsa- 
tion with period P = 202.5 ±0.5 s. Its spectrum was well fitted 
with different models (power-law, power-law plus blackbody, 
two black-bodies, cut-off power-law) and the measured flux 
implied a source luminosity of 3 x 10 34 erg s _1 between 3 and 
30 keV. RX J0440.9+4431 has been detected also in the hard 
X-ray range: it is reported (as source PBC J 0440.9+4432) in 
the P alermo Swift-BAT hard X-ray catalogue dCusumano et al 



|2010|), with a 15-150 keV flux of (2.0±1.1) x 10- n ergcm- 2 
s~\ and in the catalogue obtained with the INTEGRAL/IB IS 7- 
year All-Sky Hard X-ray Survey dKrivonos et alJEoiOa ). with 



a 17-60 keV flux of (1.36 ± 0.22) x 10~ n erg cm" 2 s" 1 . 

In the optical/IR waveband, the Ha line shows a double- 
peak profile, varying from symmetric shape to completely dis- 
torted on one side (V/R phases), with a correlation between the 
equivalent width of the Ha line and the infrared magnitudes: 
as the EW(Ha) decreases the IR magnitudes become fainter. 
This long-term optical/IR variability is attributed to structural 
changes in the Be star's circumstellar disc, which alternates de- 
cline and recover y phases on typical time scales T^ sc > 10 yrs 
dReig et alJl2005h . 

During the RossiXTE observation this pulsar showed lit- 
tle X-ray variability (by a factor < 10), therefore it was in- 
cluded in the class of persistent, low-luminosity and long- 
period Be pulsars. However, recently this source has shown 
three consecutive flux increases, spaced by ~ 5 months from 
each other: the first one was detected by MAXI/GSC on 31 



served by INTEGRAL on 1 September 2010 (IKrivonos et al 
l2010bl): the third one w as seen by Swift on 29 January 2011 



(Tsvgank ov et al.l l201 1). The mean time between the starts of 
the outbursts was ~ 155 days, and the peak luminosities of 
these flares were lower than 10 37 erg s _1 (assuming a source 
distance of 3.3 kpc). Assuming that these events occured at the 
periastron passage of the NS, where it approaches the decre- 
tion disc of the Be star, their separation is a good estimate of 
the orbital period of the system. It is interesting to note that 
the value of 155 days is also in agreement with the period of 
about 150 days der ived from the Corbet diagram of P sp i n vs 



Porut dCorbetlll986h . Therefore, after 4U 0352+308/X Persei, 
RX J0440.9+443 1 could be the second persistent Be/X-ray pul- 
sar with a known orbital period. 

2. Observations and data reduction 

RX J0440.9+4431 was observed with XMM-Newton on 2011 
March 18 (MJD = 55638 4 17). The three EPIC cameras, i.e . 
one pn dStriider et alJbOOlh and two MOS dTurneretalJl200ll) . 
were operated in Large Window mode, with a time resolution of 
48 ms for the pn camera and of 0.9 s for the two MOS cameras; 
the effective source exposure time was, respectively, of ~ 14 ks 
and ~ 17 ks. For all cameras the medium thickness filter was 
used. RX J0440.9+4431 was also observed for ~ 17 ks by the 
Reflection Grating Spectrometer (RGS), wh ich was operated in 
Spectroscopy mode dden Herder et al . 200ll). 

We used version 11.0 of the XMM-Newton Science 
Analysis System (SAS) to process the event files. After the 
standard pipeline processing, we looked for possible intervals 
of high instrumental background, with negative result. EPIC 
source events were selected within a circular area, with an ex- 
traction radius of 30" for all the cameras; the corresponding 
background events were accumulated on large circular areas 
free of sources and with radii of 120", 50" and 60" for the pn, 
MOS1, and MOS2 cameras, respectively. We selected all the 
events in the energy range 0. 15-12 keV, with pattern range 0^4 
(i.e. mono- and bi-pixel events) for the pn camera and 0-12 
(i.e. from 1 to 4 pixel events) for the two MOS. In all cases the 
background contribution to the total count rate (CR) was neg- 
ligible, resulting in a net CR of ~ 5.8 cts s _1 for the pn and ~ 
1.8 cts s _1 for each of the two MOS. Although the source CR 
was very high, we checked with the SAS task epatplot that no 
event pile-up affected our data. 

3. Timing analysis 

For the timing analysis we considered only the EPIC data. To 
measure the pulse period, we converted the event arrival times 
to the solar system barycenter and combined the three datasets. 
We me asured the pulse period by a standard phase-fitting tech- 
nique dPalFOsso et al.l 12003). obtaining a best-fit period P = 
204.96 ± 0.02 s. This result takes into account the effects of 
the source variability during the observation (Fig. [T); in fact, 
the individual pulses are detected with high significance along 
the whole observation, even when the source is faint. 
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In Fig. Q] we report the background-subtracted light curves 
in the energy ranges 0.15-3.5 (soft), 3.5-12 (hard) and 0.15- 
12 keV (total), together with the hardness-ratio HR between 
the hard (H) and soft (S) light curves (computed as H/S); the 
two energy ranges were defined in order to obtain a comparable 
number of counts, while the time bin of 205 s, corresponding 
to one pulse period, was chosen to avoid the effects due to the 
periodic pulsations. In both ranges the average CR is ~ 4.5 cts 
s~\ with an increasing trend along the observation; moreover, 
it was highly variable over the short time-scale, since there are 
CR variations up to ~ 30 % between consecutive time bins. 
Also the HR is characterized by a similar variability but, in this 
case, no long time-scale trend is observed; there is also no clear 
correlation with the source CR (Pig. O, although the fit with a 
constant value HR = 1 is rejected (x^/d.o.f. = 3.62/74). 
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Fig. 1. Background-subtracted light curves of RX 
J0440.9+4431 in the energy ranges 0.15-3.5, 3.5-12, and 
0.15-12 keV, with a time bin of 205 s (i.e. one pulse period). 
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Fig. 2. Hardness-ratio variation of RX J0440.9+443 1 as a func- 
tion of the 0.15-12 keV count rate, with a time bin of 205 s. 



The XMM-Newton observation allowed us to investigate 
the pulse profile, for the first time, even at energies below 3 



keV. In Fig.[3]we show the folded light curves in four different 
energy ranges, defined in order to obtain a comparable number 
of counts in each of them, together with the hardness ratio be- 
tween the 3.4-10 and 0.15-3.4 keV energy ranges. The shape 
of the pulse profile is similar in the four ranges: in all cases it 
shows a single broad peak, but it is not possible to fit it with a 
simple sinusoidal model; the measured pulsed fraction, defined 
as (CR max - CR min )/(2 x CR avcrag0 ), is ~ 55 % for all the 
energy ranges; the increasing part of the curve is more regu- 
lar than the decreasing one, which is steeper at the beginning 
and flatter at the end. However, we note an energy dependence 
of the pulse profile around the CR minimum (phase <fr = 0.7): 
the minimum CR value is reached at slightly later phases for 
increasing energies, and also the following CR increase is de- 
layed; moreover, at the minimum egress the CR increases more 
suddenly at high energies than at the low ones. As a conse- 
quence, the HR between the hard and soft curves is highly vari- 
able in this phase interval: it first shows a sharp minimum just 
after the CR minimum, suddenly followed by a large peak. 
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Fig. 3. Background-subtracted light curves of RX 
J0440.9+4431 in the energy ranges 0.15-2.1, 2.1-3.4, 3.4-5, 
and 5-10 keV, folded at the best-fit period P = 204.96 s. 



4. Spectral analysis 

For the EPIC source and background spectra, we adopted the 
same extraction parameters used for the light curves; we gener- 
ated the applicable response matrices and ancillary files using 
the SAS tasks rmfgen and arf gen. We considered also the 
RGS data, using the source and background spectra and the re- 
sponse matrices obtained with the standard reduction pipeline. 
In order to ensure the applicability of the x 2 statistics, the EPIC 
and RGS spectra were rebinned with a minimum of 100 and 30 
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counts per bin, respectively; then they were fitted in the energy 
range 0.3-12 keV using XSPEC 12.4.0. In the following, all 
spectral uncertainties and upper limits are given at 90 % con- 
fidence level for one interesting parameter, and we assume a 
source distance of 3.3 kpc. 

After checking that separate fits of the three EPIC cam- 
eras gave consistent results, we fitted them simultaneously to 
increase the statistics; then we checked that both RGS spec- 
tra were consistent with the EPIC ones, therefore we included 
also them in the spectral analysis; to this aim, we introduced 
relative normalization factors among the spectra of the five 
cameras. Using an absorbed power-law (PL) model, we ob- 
tained a hydrogen column density TVh = (1.19 ± 0.02) x 10 22 
cm -2 and a photon-index V = 1.16±0.01, with x^/d.o.f. = 
2.097/1083. Using an absorbed blackbody (BB) model, we ob- 
tained 7V H = (3.33 ± 0.0.9) x 10 21 cm" 2 , a BB temperature 
kT BB = 1-58 ± 0.01 keV, and a BB radius R BB = 306 ± 3 m, 
with x^/d.o.f. = 2.604/1083. 

In both cases the fits of the spectrum were unacceptable, 
with large residuals, so we repeated the fit with a PL+BB 
model. In this way we obtained a significant improvement 
(Fig. |4j, with x^/d.o.f = 1.105/1081. The corresponding best- 
fit parameters are N B = (7.3 ± 0.4) x 10 21 cm~ 2 , T = 0.85 
± 0.07, and kT BB = 1.34 ± 0.04 keV. The PL normalization 
is 7 PL = (1.8 ± 0.2) x 10~ 3 ph cm" 2 s" 1 keV" 1 at 1 keV 
and the BB radius R BB = 273 ± 16 m. The absorbed flux in 
the energy range 0.3-12 keV is / abS! x = 6.0±°; 2 x lO" 11 
erg cm~ 2 s — 1 , while the corresponding unabsorbed flux is 
/unabs.x = 6.7 x 10 -11 erg cm' 2 s~ x , which implies a source 
luminosity Lx = 8.3 x 10 34 erg s _1 . The BB component con- 
tributes for about the 35 % to the 0.3-12 keV source flux. 

Neither the EPIC nor the RGS spectra show clear emis- 
sion or absorption features. We also looked for narrow iron K a 
emission lines between 6 and 7 keV, with different widths be- 
tween and 0.5 keV. We found no evidence of such a compo- 
nent, with an upper limit on its equivalent width of ~ 70 e V (at 
90 % c.l.) at most. 

5. Phase-resolved spectroscopy 

As shown in Fig. [3] the HR is characterized by large varia- 
tions along the pulse: therefore it is interesting to investigate 
the spin-phase resolved spectrum. To this aim, we first ana- 
lyzed the background subtracted spectra of the EPIC data in 
two different phase intervals, i.e. the soft (S, <f> = 0.35-0.50) and 
the hard (H, </> = 0.75-0.90): in spite of the HR difference, they 
are characterized by a comparable CR in the total (0.15-12) 
energy range and, then, of photon counts in the accumulated 
spectra. 

The independent fits of the two spectra with a single PL 
or BB model do not provide satisfactory results (in some cases 
they are rejected by the data), while the use of a PL+BB model 
improves significantly the spectral fit goodness; the best-fit pa- 
rameters are reported in Table Q] Taking into account the es- 
timated uncertainties, none of the spectral parameters shows a 
significant variability between the S and the H spectra: in all 
cases they are characterized by comparable or consistent best- 
fit values. 
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Fig. 4. Top panel: total spectrum of RX J0440.9+443 1 with the 
best-fit PL+BB model. The spectra of the pn, MOS1 and MOS2 
cameras are shown in black, red, and green, respectively, while 
those of the RGS1 and RGS2 are shown in blue and light blue. 
Middle panels: residuals (in units of rj) between data and model 
in the case of the single PL and of the single BB. Bottom panel: 
residuals in the case of the PL+BB. 

Table 1. Best-fit spectral parameters for the phase-resolved 
spectroscopy of RX J0440.9+4431, in the case of the indepen- 
dent fit of the two spectra. 



Phase interval 


S H 


7V H (xl(F cm~ 2 ) 


0.76t u ;^ 0.69± u ;f 7 


r 


i i^+U.^b r\ QQ + U.bb 
1 - lz '-0.24 KJ - oy -0.42 

2.76+^f (63 %) 2.28+^30 (44 %) 


fcT BB (keV) 
Rbb (m) 

/bb 


1 - JJ -0.13 1 - J ^-0.11 

237± 4 7 247±p 
1.63+°-^ (37 %) 2.95+°'^ (56 %) 


/tot 
x'/d.o.f. 


4 0Q+U.18 ^1+0.22 
*- 3y -0A2 J - Z ' J -1.25 

1.15/153 1.09/170 



a Absorbed flux in the energy range 0.3-12 keV, in units of 10 1 erg 

cm -2 s _1 



We also investigated the spectral variability at the pulse 
minimum, where the HR shows an abrupt variation, by con- 
sidering the spectrum of two narrow phase ranges at <f> = 0.72 
and <j) = 0.76. Due to the limited count statistics it was not pos- 
sible to perform a detailed spectral analysis, but we checked 
that they are consistent with the best-fit models of the S and H 
phase intervals, respectively. 
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To investigate the relative variations of the two components 
with the period phase, we also simultaneously fitted the S and 
H spectra with the PL+BB model forcing the same value of 
ATh, r, and kT BB for the two phase intervals. In this case we 
obtained N B = (0.92±g-g|) xlO 22 cm" 2 , r PL = 1.48±g-& and 
kT BB = 1-67 ± 0.05 keV, with x£/d.o.f. = 1.17/328. The cor- 
responding normalization values are reported in Table|2](Case 
1). With these constraints on the model parameters, the spectral 
changes as a function of the phase are reproduced by the vari- 
ations in the relative contribution of the two components. The 
values reported in Table [2] show that, with this assumption, the 
observed spectral variability is due to both components, since 
their flux varies of more than 50 % between the two phases, 
even if in opposite directions: the PL flux decreases and the BB 
flux increases from the S to the H phase. 

In order to confirm that the thermal component varies as a 
function of the rotational phase, we should prove that a con- 
stant BB component is rejected by the data. To this aim, we 
modified the test model by linking the BB parameters for the 
two spectra together, while both the photon-index V and the 
normalization ipi, of the PL component could vary indepen- 
dently in the two phase intervals. The resulting best fit gives 
iV H = (0.84 ± 0.09) x 10 22 cm" 2 , kT BB = 1.54 ± 0.08 keV, 
and R BB = 223^^g m. The power-law parameters are shown 
in Table [2] (Case 2). Even with this model we found a good 
fit quality (x 2 /d.o.f. = 1.16/328), fully comparable to the pre- 
vious case. This result suggests that it is possible to attribute 
the whole spectral variability to the PL component, and that a 
constant BB cannot be ruled out. 



6. Discussion 

The XMM-Newton observation of RX J0440.9+443 1 is the first 
one of this source performed with a large X-ray telescope of 
the last generation; it has allowed us to investigate, for the first 
time, the timing and spectral properties of this pulsar also at 
low X-ray energies below 3 keV. 

We have obtained a new measurement of the pulse period: 
P = 204.96 ± 0.02 s. This value is in agreement with the pe- 
riod measurements obtained by the Gamma-ray Burst Monitor 
on-board the FERMI satellite since its launch in 20083 they 
show that this source is characterized by a variable spin period, 
which is typical of wind-fed binary systems. On the longer 
time-scale, our value is slightly larger than the previous value 
of 20 2.5 ± 0.5 s found in 1998 by RossiXTE dReig & Rochd 
1999b : the period difference (AP/P ~ 1%) is too large to be 
ascribed to orbital motion and implies an average spin-down 
during the past 13 years, with average P = (6.4 ±1.3) x 10 -9 
s s -1 . 



We have estimated a source luminosity Lx ~ 6 x 10 34 erg 
s" 1 in the 2-10keV energy range. To investigate the long-term 
evolution of RX J0440.9+4431, we considered also the X-ray 
flux measurements obtained with ROSAT and RossiXTE. For 
ROSATwe cons idered the CRs wh ich were obtained during the 
All Sky Survey dVoges et al. 2000h and in a subsequent, spe- 
cific observation (Motc h*et alJ 1997 ). In their case, we used the 



WebPIMMS tooE and assumed our best-fit PL+BB emission 
model to infer the source flux in the 2-10 keV energy band: 
in this way we estimated a 2-10 keV luminosity of, respec- 
tively, (1.1 ± 0.3) x 10 34 and (2.4 ± 0.2) x 10 34 erg s" 1 . 
For RossiXTE we consid ered the value obtaine d with the two- 
blackbody model used bv lReig & Rochd dl999l) . i.e. - 2 x 10 34 



erg s _1 . Although RX J0440.9+4431 was observed by XMM- 
Newton at a larger flux level, all these luminosity estimates are 
in the range 10 34 — 10 35 erg s _1 : the variation is smaller than 
a factor 10, therefore also the XMM-Newton observation sug- 
gests that RX J0440.9+4431 is a persistent BeXRB. On the 



three different events ( 


Morii et al.ll2010t Krivonos et al.ll2010b; 


Tsveankov et al. 201 1 


); based on the measured CRs, we have 



estimated that in all cases the source luminosity has risen up 
to a few 10 36 erg s . This implies a source variability of at 
least two order of magnitudes, in contrast with the classifica- 
tion of RX J0440. 9+443 1 as a persistent source. We analyzed 
also the RossiXTE ASM light-curve of RX J0440.9+443 1 since 
1996, in order to check the occurrence of other flux increases, 
but we found no clear evidence of such type of events in addi- 
tion to the previous three ones; therefore the flux level of RX 
J0440. 9+4431 is usually consistent with a persistent nature of 
the source, while the transient behaviour recently observed is 
likely due to structural changes in the circumstellar disc of the 
Be star. 

From the spectral analysis, we obtained a hydrogen col- 
umn density N B = (7.3 ± 0.4) xlO 21 cm -2 , which is much 
lower (almost one order of magnitude) than the values esti- 
mated by RossiXTE. However, we note that the energy range 
of the RossiXTE spectral analysis (above 3 keV) is not well- 
suited for a good estimate of A^h, while the low-energy end 
of the XMM-Newton spectra (0.3 keV) allows a much more 
reliable analysis. LS V +44 17, the optical count erpart of RX 
J0440.9+4431, has a color excess E(B-V) = 0.65 dReigll201 ll) . 
Assuming Ay = 3.1 E(B-V) and the average relation Ay = 
N B x 5.59 x 10~ 2 2 cm" 2 between opti cal extinction and 



X-ray absorption (Pred ehl & Schmittlll995l) . this would give 
iVn = 3.6 x 10 21 cm" 2 , a value which is a factor ~ 2 lower than 
our result. We note that the data used to calibrate the previous 
relation is characterized by a comparable dispersion, therefore 
there is a rough agreement between the two measurements. 

We found evidence of a previously undetected thermal 
component, in addition to the main power-law: it is character- 
ized by a high temperature (fcT = 1.34 keV) and a small emis- 
sion area (R = 273 m), and contributes for ~ 35 % of the source 
flux below 10 keV. On the other hand, it was not possible to fit 
the XMM-Newton spectr a with the two-blackbody model used 
by iReig & Rochd (11999b . The RossiXTE spectrum showed a 
possible feature at ~ 6.2 keV, suggesting the presence of a Fe- 
K emission line; but it was not possible to constrain its param- 
eters, and only an upper limit of ~ 100 eV on the Equivalent 
Width (EQW) could be provided. In spite of the higher lumi- 
nosity level, we found no evidence of the Fe-K emission line 
and we set an upper limit of ~ 70 eV (at 90 % c.l.) to its EQW. 



1 http://gammaray.nsstc.nasa.gov/gbm/science/pulsars/ligfa^ 
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Table 2. Best-fit values for the black-body and power-law parameters, when the two spectra are fitted simultaneously with 
common values of JVh, I\ and /cTbb {case 1) or with common values of JVh, kT&&, and Rbb {case 2). 



Phase 
interval 


case 1 
S H 


case 2 
S H 




n n ? +U.U8 n Q,+u.ua 

"• yz -0.07 "- yz -0.07 


0.84 ± 0.09 0.84 ± 0.09 


r 

T b 

J PL 
fpL 


148 +>i - 17 148 M>LV 

o o 7 +0.55 9 9/1+0.49 
3'3 l _ q_46 ^"^+_041 

2.34±°;g (54 %) 1.56±g;gi (31 %) 


1 - JZ -0.19 1 - uu -0.15 
9 /r S +0.68 , no+0.44 
—0.59 —0.40 

1.73lg;^ (42 %) 2.93±S:2e < 55 % ) 


kTsB 
Rbb 
/bb 


i.67± 0.05 1.67 ±0.05 
2.02 +0 % (46 %) 3.47+2'oe (69 %) 


J.54 ± 0.0S 7.54 ± 0.08 

990+17 993+17 
ZZJ -16 ZZJ -16 

2.42+°'°« (58 %) 2.42+°-^ (45 %) 


/tot 

xl 

d.o.f. 


4 36 +uilJ 5 03 +ui4 

H - ,o -0.16 J - u:) -0.11 

1.17 

328 


4,^+u.ilj s «+u.uu 

H - 1J -0.20 J "-> J -0.28 
1.16 

328 



a xlO 22 cm" 2 

6 x 10~ 3 ph cm" 2 s" 1 keV" 1 at 1 keV 

c Absorbed flux in the energy range 0.3-12 keV, in units of 10 -11 erg cm -2 s 



For its size and temperature the BB excess observed in 
RX J0440.9+4431 is similar to those detected in the other 
persistent Be/NS XBPs, i.e. LS I +6 1° 235/RX J0146.9+6121 



(|La Palombara & Mereghettil l2006h. X Persei/4U 03 52+309 
dCoburn et all l2Q0lt lLa Pa lombara & Mereg hettil 
and LS 1698/RX J1037.5-5647 dReig & Rochd 



2007) £ 



1999, 3 



La P alombara et al. I l2009h : therefore this hot BB spectral com 
ponent is a common property of this type of sources. However, 

we note that a similar feature has been observed also in 3A 

1 11 1 1 11 

0535+262 (Orlandini et al. 2004; Mukheriee & Paul 2005), 4U 



2206+ 54 dMasetti et all 12004: Torreion et al. 2004; Re ig et al 
2009|), and SAX J2103.5+4545 dlnam et alj|2004l) which are 



other types of XBPs with low luminosity (Lx < 10 erg 
s" 1 ) and long pulse period {P > 100 s). Moreover, recently 
the same type of BB excess has been detected also in the three 
Supergiant Fast X -ray Transients (SF XTs) IGR Jl 1215- 5292 
dSidoh et alJl2007j). IGR J0840 8-4503 dSidoh et alJl2009h and 
XTE J1739-302 jBozzo et alj l2O10b . only the first of which 
is a confirmed pulsar {P = 187 s). Finally, very recently a 
hot BB excess has been detected also in the SMC binary 
pulsar SXP 1062, a possible new persistent Be X-ray binary 
dHenault-Brunet et al. 201 lh . In Fig. [5] we report the best- 
fit radius and temperature for the BB component of these 
sources, together with lines showing four different levels of 
the blackbody luminosity; for some source more than one set 
of values is shown, corresponding to different observations or 
flux levels. In most cases the spectral parameters are within 
a narrow range of values, i.e. A;Tbb ~ 1-2 keV and i?BB < 
200 m: we emphasize that, in all these cases, the estimated 
total source X-ray luminosity is ~ 10 34 erg s _1 , with a 
20^40 % contribution of the blackbody component. For RX 
J0440.9+4431 the BB radius is slightly higher (Rbb ~ 270 m), 
in agreement with the fact that its total X-ray luminosity is also 
higher [Lx ~ 8 x 10 34 erg s" 1 ): when observed at high (i.e. 
L x ~ 10 35 erg s" 1 ) luminosity level, also RX J1037.5-5647 
(point 2), 4U 0352+309 (point 4), SAX J2 103 .5+4545 (point 
10) and IGR Jl 1215-5292 (point 11) show large values of 
temperature and/or radius. 




200 500 
Blackbody radius (m) 

Fig. 5. Best-fit values for radius and temperature of the 
BB component in the case of RX J0440.9+4431 (empty 
star), RX J1037.5-5647 (filled circles), 4U 0352+309 (filled 
squares), RX J0146.9+6121 (empty circle), 3A 0535+262 
(empty square), SAX J2 103. 5+4545 (empty triangle), 4U 
2206+54 (filled triangles), IGR Jl 1215-5292 (asterisks), IGR 
J08408-4503 (cross), XTE J 1739-302 (filled star), and SXP 
1062 (empty diamonds). The continuous lines connect the 
blackbody parameters corresponding to four different lev- 
els of_Jurmnosity_(in erg_s_~ 1 ). References: - this w ork; 
1-lLa Palombara et alj d2009b: 2 - iReig & Rochd d 1999b: 3 - 
ICoburn et al.1 d200lb: 4 - lLa Palombara & Mereghettil d200 / 
5 - [Mukher e e & Paull d2005b: 6 - lLa Palombara & Mereghett 
Masetti et alj d2004: 8~lTorrei6n et al.1 d20Q. " 



d2006h: 



7 - 



Reig et al. 



( 2009): 10 - linam etalJ J2004l>: 11 



iHenault-Brunet et alJ d201 lb" 



J2007b: 12 -ISidoli et alJ d2009b : 13 -Ibozzo et alJ d2oTob: 14 - 



Sidoli et al 



In contrast to this sample of sources, several XBPs are char- 
acterized by a soft excess, since the fit of this component with a 
thermal emission model provides low temperatures (kT < 0.5 
keV) and large emitting regions (R > 100 km). In Fig. [6] we 
report the luminosity and pulse period of both types of XBPs: 
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the soft excess and the hot BB ones are reported as squares and 
circles, respectively. Based on their distribution in the P — Lx 
diagram, these pulsars are divided into two distinct groups: the 
sources in the first group are characterized by high luminos- 
ity (L x > 10 37 erg s _1 ) and short pulse period (P < 100 s), 
and in most cases they are in close binary systems with an ac- 
cretion disk; those in the second group have low luminosities 
(Lx < 10 36 erg s -1 ) and long pulse periods (P > 100 s), since 
they have wide orbits and are wind-fed systems. While all the 
pulsars in the first group are characterized by a soft excess, both 
types of pulsars are present in the second group. In this case, 
the hot BB pulsars are the ones that, on the average, are charac- 
terized by the lowest luminosities and the longest periods. This 
suggests that the hot BB spectral component is a common fea- 
ture of the low-luminosity and long-period XBPs. However, 
in the second group of sources there is no clear separation be- 
tween the two types of pulsars, since there is a partial overlap 
between the soft excess and the hot BB ones. Therefore, the ori- 
gin of the spectral difference in pulsars having comparable val- 
ues of luminosity and pulse period is unclear. Even if it could 
be explained by really different emission mechanisms at work 
in these sources, we can not exclude that, at least in some cases, 
it is due to the degeneracy of spectral fitting in some spectra. In 
fact, on the basis of the X-ray spectrum it is often impossible 
to distinguish between a BB with small radius and high temper- 
ature and one with large radius and low temperature (see, e.g., 
Bozzo et al.ll2oTob . 



□ □ 



cP 



33 
0.01 



1 10 

Pulse Period (s) 



Fig. 6. X-ray luminosity (in the 2-10 keV energy range) of the 
XBPs with a detected thermal excess as a function of the pulse 
period. The filled circle refers to the XMM-Newton observation 
of RX J0440.9+4431, empty circles refer to other detections of 
hot BB pulsars, open squares to the soft excess sources. 



Based on our result s and on the emission models proposed 
bv lHickoxetail d2004 . also in the case of RX J0440.9+4431 
it makes sense to attribute the observed spectral excess to the 
thermal emission from the NS polar cap, in agreement to what 
already suggested for the other persistent Be pulsars. In or- 
der to check this hypothesis, we assume that the source is in 
the 'accretor' status, with matter accretion on the NS surface. 



= 10 cm, the source lu- 
1 implies an accretion rate 
10 12 G, 



Assuming Mns =1.4 M© and i?Ns 
minosity Lx = 8.3 x 10 34 erg s 

M~4.5xl0 14 gs _1 and, adoptin g 7j NS = 10 12 G, a m ag- 
netospheric radius R rn ~ 10 9 cm ( Campana et alJ fl998l) . In 
this case, based on the relation R co i ~ i?NS (RNs/Rm) ' 5 
dHickox et al. 2004), we would obtain R co i ~ 320 m. This 
value, which is an estimate of the expected size of the polar 
cap, is remarkably near to the estimated blackbody emitting 
radius (Rbb ~ 270 m): therefore it strongly supports the sug- 
gestion of a polar-cap origin for the observed BB emission. If 
this description is correct, we would expect to observe some 
variability of the thermal component along the pulse phase. 
The phase-resolved spectral analysis confirmed the spectral 
variability along the pulse, but it can be attributed to the PL 
component, since a constant BB component is fully compat- 
ible with the spectral data. The lack of variability of the BB 
component is not in contradiction with a polar-cap origin, as 
an anisotropic radiation can appear steady along the pulse for 
some ge ometrical configura tions of the hot spots on the NS 
surface dBeloborodovtf2 002). On the other hand, the phase de- 
pendence of the non-thermal component is likely due to the 
anisotropy of the Comptonized radiation of the accretion col- 
umn . This is shown , for ex ample, in the self-consistent model 
bv lBecker& Wolfjd2007h for the dynamics and the radiative 
transfer occurring in the accretion column of bright X-ray pul- 
sars. This model is based o n the physical picture originally pro- 
posed by iDavidsonl d 1973b . in which the accreting gas passes 
through a radiative, radiation-dominated shock before settling 
onto the NS surface. For the relatively low luminosity of RX 
J0440. 9+443 1 the shock could be absent or have only a limited 
extent above the neutron star surface. Thus the blackbody emis- 
sion from the dense thermal mound at the base of the column 
could be the dominant seed component for the Comptonization. 
In this geometry, the radiation is probably emerging as a 'pen- 
cil beam' rather than from the column's lateral surface as in the 
Becker & Wolff model. 



7. Conclusions 

We have analysed a ~ 17 ks XMM-Newton observation of the 
Be/NS X-ray pulsar RX J0440.9+443 1 . The source was de- 
tected at a luminosity level Lx — 8.3 x 10 34 erg s _1 in the 
0.3-12 keV energy range: this value is less than one order of 
magnitude higher than the average level of the previous ROSAT 
and RossiXTE observations, therefore it confirms the persistent 
nature of this source. This classification is not denied by the 
large flux increases recently observed, which imply luminos- 
ity levels of a few 10 36 erg s _1 , since they can be explained 
by structural changes in the circumstellar disc of the Be mass 
donor. 

The XMM-Newton observation has provided a refined 
pulse period P = 204.96 ± 0.02 s, which, compared to previ- 
ous measurements, implies an average pulsar spin-down P = 



(6.4 ± 1.3) x 10~ 9 s s" 1 between 1998 and 2011. It implies a 
low average accretion rate, in agreement with the persistent na- 
ture of the pulsar. The pulse profile shows a complex structure, 
which is not sinusoidal, and at all energies the PF is ~ 55 %. 
However the HR is characterized by significant variations, ex- 
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pecially around the pulse minimum, which indicates a spectral 
variability of the source. 

The source spectrum shows a count excess above the main 
power-law component, described by a blackbody with kTsB — 
1 .3 keV and R B b — 270 m, which contributes ~ 35 % to the 
source luminosity between 0.3 and 12 keV. We found no evi- 
dence of a narrow iron K Q line between 6 and 7 keV, with an 
upper limit of ~ 70 eV on its equivalent width. The blackbody 
radius is comparable to the estimated size of the NS polar-cap, 
which suggests that the origin of this component is on the NS 
surface. The phase-resolved spectroscopy neither confirms nor 
disproves this scenario, since the observed spectral variability 
along the pulse period can be attributed to the power-law com- 
ponent and is consistent with a constant thermal component. 

The spectral properties of RX J0440.9+4431 are in full 
agreement with those observed in the other three persistent 
Be binary pulsars 4U 0352+309, RX J0146.9+6121 and RX 
J1037. 5-5647, thus confirming that the hotBB spectral compo- 
nent is a common property of this type of sources. Moreover, 
we have shown that the same type of feature has been detected 
also in other low-luminosity and long-period pulsars; therefore 
it is an ubiquitous phenomenon which requires further investi- 
gations. 
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